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Abstract 

By choosing a dynamical system with d different couplings, one can rearrange a system 
based on the graph with given vertex dependent on the dynamical system elements. The 
relation between the dynamical elements (coupling) is replaced by a relation between the 
vertexes. Based on the Eq transverse projection operator. We addressed synchronization 
problem of an array of the linearly coupled map lattices of identical discrete time systems. 
The synchronization rate is determined by the second largest eigenvalue of the transition 
probability matrix. Algebraic properties of the Bose-Mesner algebra with an associated 
scheme with definite spectrum has been used in order to study the stability of the coupled 
map lattice. Associated schemes play a key role and may lead to analytical methods in 
studying the stability of the dynamical systems. The relation between the coupling 
parameters and the chaotic region is presented. It is shown that the feasible region is 
analytically determined by the number of couplings (i.e, by increasing the number of 
coupled maps the feasible region is restricted). It is very easy to apply our criteria to the 
system being studied and they encompass a wide range of coupling schemes including 
most of the popularly used ones in the literature. 

Keywords: Synchronization, coupled maps, dynamical systems and Bose-Mesner 
algebra. 

PACS index. 05.45.Xt,84.35.Xt 
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1 Introduction 

Simulation of the natural phenomena is one of the most important research fields and coupled 
map lattices are a paradigm for studying fundamental questions in spatially extended dynami- 
cal systems. The early definition of the coupled map lattice goes back to the Kaneko's paper [1]. 
Now we can divide them in two categories: internal and external coupled map lattice [2], Glob- 
ally coupled map is one of the most well know examples of the external coupling, with different 
number of elements[3]. Variety of couplings such as the week coupling[4], noisy coupling[5] and 
functional coupling[6] used to made the new coupled map. 

On the other hand many areas of research lie on the interface of algebra and physics. Bose- 
Mesner algebra arose independently in three areas: statistical designs, centralizer algebras of 
permutation groups, distance-transitive graphs and quantum walk [7, 8, 9, 10, 11]. Associ- 
ated schemes are combinatorial objects that allow us solve problems in several branches of 
mathematics. They have been used in the study of permutation groups and graphs and also 
in the design of experiments, coding theory, partition designs etc. In this paper the authors 
are interested to use the ability of algebraic properties of the Bose-Mesner algebra associated 
with an association scheme with known spectrum in order to study the stability of the coupled 
map lattice, since elements of Bose-Mesner algebra have capability that can be diagonalized 
simultaneously. 

The aim of the present paper is threefold: introducing a generalized model of the coupled 
map lattice which covers internal and external couplings in a form of associated schemes. By 
taking into account that idempotents are another basis of this algebra (E a longitude projec- 
tion operator, Ep transverse projection operator), in our proposed model all of the elements 
have a common prolongation. Idempotents are independent of coupling and coupling method 
influences the stability of the coupled map lattices. 

We addressed synchronization problem of an array of the linearly coupled map lattices of 
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identical discrete-time systems. We present the proof that the synchronization rate is mainly 
determined by the second largest eigenvalue of the Floquet multipliers. We also prove that 
topological behavior can be studied by Sinai-Bowen-Rule (SBR) measure for the coupled map 
lattice of measurable dynamical system at the complete synchronized state. 
In order to detect the feasible region, the stability analysis of the coupled map lattice with 
regard to the Lyapunov exponent is possible. We define the Lyapunov function by refereing to 
the projection operators E^s in order to confirm the stability analysis based on the Lyapunov 
exponents. Lyapunov function takes value at synchronized state. We have investigated some 
conditions on transition from regularity to chaos with Li-yorke theorem in a coupled map net- 
work, whose individual nodes are non-chaotic before being connected to the network. It has 
been discovered that, such a transition to chaos depends not only on the network topology but 
also on the original node dynamics. 

In order to determine whether our results in this paper are generic or they are only limited to 
the specific model, we study several other coupled maps. It was proved that the global coupled 
map lattice corresponds to the compete graph; and the nearest neighbor coupled map can be 
studied based on the model of Associated schemes. It is very easy to apply our criteria to the 
system being studied, and they encompass a wide range of coupling schemes including most 
of the popularly used ones in the literature. 

The paper is organized as follows: In section 2 of the paper we present a brief outline of some 
of the main features of the associated schemes, such as adjacency matrix, distance- regular 
graphs, stratification and orthonormal basis of the strata. In Section 3, the generic model 
of the discrete-time coupled map Lattice (Internal and coupled map lattice) is introduced in 
adjacency matric perspective. In section 3, the synchronization conditions are also studied. 
Section 4, discusses the existence of SBR-measure at the synchronized state in the light of 
a measurable dynamical system. In the section 5, the stability analysis of the coupled map 
lattice based on the discussion of Floquet multipliers of coupled map lattice is done through 
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three subsections; Lyapunov exponent calculation, Lyaunov function analysis and Li-yorke 
theorem. In section 6, we present some examples of the coupled map lattice model and their 
corespondent graph such as the complete graph, strongly regular graphs and cycle graph. This 
section followed by an outlook section. Respective appendix can be found at the end of paper. 

2 Association schemes 

In this section we give a brief outline of some of the main features of association scheme, and the 
reader is referred to [7] for further information on the Associated schemes. Recall that a finite 
graph is a finite set I, whose elements are called vertices described by vertex set V, together 
with a set of 2 — subset of T called edges. An association scheme with d associate classes on a 
finite set T is coloring of all the edges of the graph by s colors such that the set of edges with 
the same color form a non empty relation subsets on T (i,e, subset of T x F ), denoted by R iy 
where i — 1, ..,d corresponds to different colors. Then the pair Y = (T, {Ri}o<i<d) consisting 
of a set T and a set of relations {Ri}o<i<d together with the following 4 conditions is called an 
association scheme. 

1- {Ri}o<i<d is a partition of T x T 

2- Rq = {(a, a) : a G T} 

3- Ri = R\ for < % < d, where R\ = {(/?, a) : {a, (3) e R { } 

4- Given (a,P) £ R k ,Pij =| {7 G T : (a, 7) G Ri and (7,/?) G Rj} |, 

where the constants are called the intersection numbers, depending only on i,j and k 
and not on the choice of (a, (3) G Rk- Then the number N of the vertices V is called the 
order of the association scheme and Ri is called a relation or associate class. Let T = (V, R) 
denote a finite, undirected, connected graph, with vertex set V edge set R path-length distance 
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function d, and diameter d := max{d(a, (3) : a, (3 e V}. For all a, (3 G V and all integer i ,we 
set Ft = {V, Ri) = {(a,P):a,/3 eV : d(a, (3) = i] so that T^a) = {f3 e V : d(a, (3) = i}. 

3 The Bose-Mesner algebra 

Let R denote the field of complex numbers. By Matv(R) we mean the i?-algebra consisting of 

all matrices whose entries are in R and whose rows and columns are indexed by V. For each 

integer i (0 < a < d), let A a denote the matrix in Maty{R) with (i,j)-entry 

f 1 if (i,j) E R a , 

(A a )i, 3 = (iJeT). (3.1) 
1 otherwise 

The matrix A a is called an adjacency matrix of the association scheme [8]. Then we have 
A = /(by (2) above) and 

A a A p = itpl^, (3.2) 

7=0 

(by (4) above), so A Q , Ai, A d form a basis for a commutative algebra A of Mat v (C), where 
A is known as the Bose-Mesner algebra of Y . Since the matrices commute, they can be 
diagonalized simultaneously, that is, there exists a matrix S such that for each A e A, S _1 AS 
is a diagonal matrix. Therefore A is semi-simple and has a second basis E ,...,E d (see [9]). 
These are matrices satisfying 

1 d 

Eq = JjJN, E a Ef3 = 5 a f3E a , ^ E a = I N . (3.3) 

The matrix j^Jn (where J N is the all-one matrix in A and iV =| T \ ) is a minimal idempotent 
(idempotent is clear, and minimal follows since rank matrix Jn equal 1. The E a , for (0 < 
a,f3<d) are known as the primitive idempotent of Y. Let P and Q be the matrices relating 
our two bases for A: 

d 

A P = J2 P ^ 0<P<d, (3.4) 

1 d 

E P = Tf E Q«pA a , 0<f3<d. (3.5) 
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Then clearly 

PQ = QP = NI N . (3.6) 

It also follows that 

ApE a = PapEa, (3.7) 

which shows that the P a p (resp. Q a p) is the cn-th eigenvalues (resp. the a-th dual eigenvalues 
) of Ap (resp. Ep) and that the columns of E a are the corresponding eigenvectors. Thus 
m a = rank(£' Q ,) is the multiplicity of the eigenvalue P a p of Ap (provided that P a p ^ P 7j a for 
7 7^ a). We see that m = 1, J2 a m a = N, and m a =traceE a = N(E a )p/3 (indeed, E a has only 
eigenvalues and 1, so rank(i? 7 ) equals the sum of the eigenvalues). Also, the eigenvalues and 
dual eigenvalues satisfy 

rripPp a = k a Q afS , 0<a,(3<d, (3.8) 

where for all integer a (0 < a < d), set k a = p° aa = (N — l) a C%, and note that k a ^ 0, since 
R a is non-empty. 

3.1 Definition of generic N-coupled map 

Coupled map lattices are arrays of states whose values are continuous, usually within the unit 
interval, or discrete space and time. The generic N-coupled map dynamical system can be 
considered as a N-dimensional dynamical map defined as: 



$(xi(ri), ...,x N (n)) 



x 1 (n+ 1) = Fi(xi(n), ...,x N (n)) 
x 2 {n + 1) = F 2 (x 1 (n), ...,x N (n)) 



(3.9) 



^ x N (n + 1) = F N {x 1 (n), ...,x N (n)). 

Adjacency matrix element such as A can cover the elements of the dynamical systems in the 
coupled map lattice. Ai, A 2 , etc, have the capability of generating the different coupled map 
matrices such as nearest neighbor coupled map, second nearest neighbor coupled map and 
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globally coupled maps. Definition of the coupled map model based on the adjacency matrix 
allows us to use the spectrum of graph based on the associated scheme. Therefore, generic 
N-coupled map is defined as follows: 



Xi(n + 1) 



a=0 K " Vxe^fxifnn / \ Q=0 z6r a (xi(n)) 



. (3.10) 



^x£r a (xj(n)) 

Where n presents the time, N is the number of the coupled map, e corresponds to internal and 
external-coupling, e^(e^) stand for coupling parameters and $ is a ID map in this model. For 
e ' a ( e a) ~^ 0; there is no coupling at all, hence, local neighborhoods have no influence on the 
behavior of the coupled map lattices. Using Bose-Mesner algebra language, we can write: 



xAn + 1) 



where 



«E4*W»)) + 1 1 - E )*E4>») 



a oc a "'a 



, i,j = l,..,N 



(3.11) 



(3.12) 



A' and A" are the elements of Bose-Mesner Algebra. e' a and e a with the same relation in the 
associated schemes are the coupling constant in the coupled map lattice 

Ei = i. E£ = i (3-i3) 

a=0 a=0 

Dynamical system symmetry consists of permutations which do not affect the adjacency ma- 
trices and elements of Bose-mesner algebra. 



a- By considering e = 0, we can generate Internal-coupled maps as follows: 



x 



(n+i) = *(E4*» = * Er( E *(*)) 



(3.14) 



b- By considering e = 1, we can generate External-coupled maps as follows: 



Xz (n + 1) = e 4 $ (**(»)) = Ef E 



(3.15) 
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3.2 Local and globally coupled maps 



We generalize the models (3.11) to the Coupled map lattices with local-global couplings on a 
one dimensional lattice of length N with periodic boundary conditions[14]. These maps follow 
Eq. (3.11) if we modify the coupling constants in the following way: 



a = 0, 1, d 



k a k a N 

where e a and 7 are coupling parameters. By substituting e a in (3.11): 



(3.16) 



x 



V a=l I a=l \xer a (xi(n)) , 



+(l-e)$ 



;i-EC-7)^)+(l-^)Er( E *(«)) 

a=l a=l xer„(xi(n)) 



+ 



_7 
TV 



«E( E ^(x(n))) + (l-c)$(X:( E *("))) 

a=0 xer a (xi(n)) a=0 xer a (x 4 (n)) 



(3.17) 



3.3 Synchronization of coupled maps 

Synchronization of two (or more) chaotic dynamical systems (starting with different initial 
conditions) means that their chaotic trajectories remain in step with each other during the 
temporal evolution. Complete synchronization was popularized after the seminal papers of 
Pecora and Carroll [12, 13]. 

Complete synchronization in the coupled map means the existence of an invariant one-dimensional 
sub-manifold x\ = ... = xn or (xi(n) = ... = xn(ji) x\{n + 1) = ... = x^in + 1)). 
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As it was discussed in section 3, the matrices A' and A" appearing in coupled map (3.11) 



are semi-simple, and hence possess the minimal idempotent E , E d , where 





) 



(3.18) 




projects coupled dynamical system on synchronized state. The remaining ones, on the other 
hand, will project it on the transverse modes. By choosing i — j and considering (3.12) the 
generic model of coupled map lattice Eq. (3.11) is reduced to Xi(n + 1) = $(xj(n)). Thus, 
at synchronized state, the behavior of coupled map lattice is studied by the behavior of single 
dynamical element. Here in this article we study the synchronization of the coupled map 
lattices with different coupling topologies such as the global coupling, the nearest neighbor 
coupling and the coupling varying with distance etc. 

In this paper we are interested in explaining the generic model based on the hierarchy of one 
and many-parameter chaotic maps which are introduced in our previous papers[17, 18, 19]. 
One-parameter families of chaotic maps of the interval [0, 1] with an invariant measure can be 
defined as the ratio of polynomials of degree N (See [15] for more detail): 



If F is substituted with chebyshev polynomial of type one T N (x), we will get & N '(x,a) and 
if F is substituted with chebyshev polynomial of type two Un(x), we will get ^\x,a). In 
hierarchy of elliptic chaotic maps F is substituted by Jacobian elliptic functions of cn and sn 
types. Also, we present some examples of discrete-time one dimensional dynamical system in 
Appendix B. 




(3.19) 
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4 Invariant measure at synchronized state 

The measure which describes the ergodic properties with respect to the typical initial conditions 
is usually called SRB measure [20]. The difficulty in proving rigorously that a given coupled 
map lattice exhibits spatio-temporal chaos lies in the finding of such a SRB measure, which 
has the following important properties: 

a- The measure is invariant with respect to symmetry transformation of N-dimensional dy- 
namical system, namely it corresponds to the scalar representation of its symmetry group. 

b- The measure is smooth along unstable eigendirections in the phase space. 

c- The measure has strong ergodic properties, including mixing and positive KS-entropy[21]. 

Symmetric transformation are kind of transformations which leave the adjacency matrices 
invariant. As we will see later, by choosing the one-dimensional maps with an invariant 
measure, such as logistics map or the ones introduced by authors in [17, 18], the coupled 
maps can display the invariant measure at synchronized state. The corresponding Forbenious- 
Perron (FP) equation for N-coupled map (3.9) is shown: 

fi(x 1 (n+ 1), ...,x N (n+ 1)) = J dx 1 ... J dx N 5 (x 1 (n+ 1) - Fi(x 1 (n), ...,x N (ri))) ... 

5(x N (n + 1) - F N (xi(n), x N (n)))^(xi(n), .., x N (n)) (4.1) 

Considering the above mentioned properties of the invariant measure, it is obvious that, at 
synchronized state it should have the following form: 

/i(xi, ...,x N ) = 5{x 2 - Xi)....6(x N - xi)/i(a;i), (4.2) 

where fi(xi) corresponds to the invariant measure of one-dimensional map defined in one- 
dimensional invariant manifold x\ — x-i — ... — xn- By exercising the synchronization condi- 
tion in Eq. (4.1), we have: 

fi(xi(n + 1), ...,x N (n + 1)) = / dxi... / dx N 5(xi(n + 1) - Fi(xi(n), x N (n)))... 
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5(x N (n + 1) - F N (xi(n), x N (n)))S(x 2 (n) - xi(n))....S(x N (n) - xi(n))^(xi) 

= J dXl MMn + 1) - *(*,(.,), ..,*,(»)))...*(**(» + D " ^(*.(»), *.(»))M*d 
= 5(x 2 (n+l)— xi(n+l))..5(xjv(n+l)— xi(n+l)) J dxiS(x 1 (n+l)—F 1 (xi(n) : .., xi(n)))/i(xi(n)), 

Now, if the one-dimensional map x(n+ 1) = xi(n)) possesses the invariant measure 

/i(xi(n)), we can write: 

H(x(n + 1)) = J 5(x(n + 1) - a^ra))^^), (4.3) 

Then, we have 

Hix^n+l), ...,x N (n+l)) = 5(x 1 (n+l)-x 2 (n+l))....5(x N (n+l)-x 1 (n+l))fj,(x 1 (n+l)), (4.4) 

We have already derived analytically invariant measure for One-parameter families of chaotic 
maps (3.19) by using arbitrary values of the control parameter a and for each integer values 
of N. 

Jx,P) = - ; — , (3>0 (4.5) 

^■W'^ vr y^IT^) (/3 + (i _ ^ " 1 ; 

It should be mentioned that reader may refer to our previous papers for derivation and the 

relation between the control parameter and (3 [18, 19, 20]. 



5 Stability analysis at synchronized state 

In this section we present the stability analysis of the coupled map lattice synchronized state 
of fix point type (i.e, one-dimensional map (3.19) at synchronized state posses fix point) and 
synchronized chaotic type ( one-dimensional map at synchronized state is chaotic and posses 
an invariant measure). This section also studies Li-Yorke type of chaos as one-synchronized 
state becomes unstable. 
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5.1 Stability analysis by Lyapunov exponent spectra 

The stability of the coupled map can be assessed by computing its Lyapunov exponent spec- 
trum. The spectrum of Lyapunov exponents of coupled map lattice with respect to the syn- 
chronization state can be evaluated in a way similar to that of one dimensional local maps 
[24]. At synchronized state, the Lyapunov exponents Ap of N-dimensional dynamical system 
described by the map (3.11) can be obtained by perturbing the state x,(n + 1) infinitesimally 
around the synchronized state (xi(0) = x, ...,x = xjv(O)); one can show that 5x(m + 1), the 
corresponding perturbed state at the time m, becomes [20]: 

5xi{m + 1) = XX ~7T~7 + ^ )x 1 =...=x N Sx i {m), (5.1) 
■ OXj [m) 

Then, by taking derivative of Xi(m + 1) with respect to Xj(m), in equation of (3.11), and 
exercising the result in (5.1), we obtain: 

5x(m + 1) = &{x{m)){eA' + (1 - e)A")8x(m) (5.2) 

by iterating (5.2) we obtain: 

n-1 

5x(n) = ( J] &(x(m)))(eA' + (1 - e)A") m 6x{0), (5.3) 

m=0 

with 

V = E fPp* (5-4) 

where r] a = ee' a + (1 — e)e' a . By using equation (3.5) and idempotency property of E a , we 
obtain: 

substituting in (5.3): 

5x(n) = n 6(x(m)) ^ \%Ep6x(0), 

m=0 (3 
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where 

rt— 1 d n—l 

6x{n) = J] 4>(x(m))E 8x(0) + ®(x(m))\™E p 5x(0) (5.5) 

m=0 0=1 m=0 

It should be reminded that in Bose-mesner algebra perspective, by resorting to projection 
operator, one could diagonalize the adjacency matrix, where diagonalization is the common 
problem in the coupled map lattice. 

E 8x(0) is a prolongation which depends on synchronization state whereas Ep8x(0) (/? = 
l,...,d) prolongations depend on transverse modes. The eigenvalues (ni^o^^N)^™) are 
called the Floquet (stability) multipliers of the orbit. The dynamical system phase space is 
formed by sum-direct of sub-space due to its idempotents project. Thus, for stabilities of 
subspace, the absolute value of Floquet multipliers should be less than 1. 

n—l n—l 

po = I] k*(m)) P p=U *(x(m))X? ((3 = 1,.-., d) (5.6) 

m=0 m=0 

Note that po is eigenvalue which corresponds to synchronous periodic orbit while pp to trans- 
verse states. The eigenvalues of the Floquet (stability) multipliers are real having only 1 value 
in magnitude . We will show them in non-increasing order: 

I PO |>| P2 |> ••• >| Pd | 

As the Eq. (5.4,6) shows, pp directly corresponds to Ep (projection operator). So, by con- 
sidering the Eq. (3.6), one can diagonalize adjacency matrix spontaneously. This simplifies 
the calculation of pp. Stability analysis of the coupled map lattice based on the stability of 
synchronous periodic orbit and transverse state, should be done in one of these three categories: 

a- Both of synchronous periodic orbit and transverse states are stable. So, feasible region is 
studied by Lyapunov exponent and Lyapunov function. 

b- Synchronous periodic orbit is unstable, but transverse states are stable. So, feasible region 
is studied by Lyapunov exponent and Lyapunov function. 
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c- Both synchronous periodic orbit and transverse states are unstable. So, stability analysis 
is studied by Li-yorke theorem. 

By taking the logarithm of power ^ absolute value of Floquet multipliers, one can study analyt- 
ically the stability of dynamical systems with Lyapunov exponent, whereas Floquet multipliers 
should be studied numerically. Notice that, 5x(0) belongs to corresponding subspace of E a , 
So, the Lyapunov exponents of N-dimensional dynamical system are defined as: 

A«= lim -]n( l^^l l =ln| A* I + lim - V In \&(x(m))\ (5.7) 

By considering one-dimensional ergodic map as a dynamical element, the last part of the 
above-mentioned Eq. (5.7) is reduced to: 

lim - Y, ln \®'( x ( m ))\ = / d(n(x)) In \&(x(m))\ = A L ($) (5.8) 

n ^°° 71 rn=0 J° 

where A^ shows Lyapunov exponent of one dimensional map. Therefore, the Lyapunov expo- 
nents of N-dimensional dynamical system is defined by: 

A /3 = ln | \p | +A L ($). (5.9) 

If in the Eq. (5.9) is substituted by "0", the equation would be reduced to: 

A o = Er p o a = i (s.io) 

hence, A = Al(<&). According to item b, for stability of transverse modes it is necessary to 
Ap < 0, (P — 1, ...,d) therefore: 

IE^I<^ AlW (5.H) 

a en 

as it was discussed Ys^oVa = 1- Therefore, synchronized state makes the coupling constants 
meet the following inequality, condition: 

1 _ e -M*) < J2(k a - < 1 + e~ x ^ {3 = 1, .., d (5.12) 

a=l ^ a 
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By considering that the multiples of rj a are positive, the gradient of this plans is in the first 
region. So, the Eq. (5.12) defines a hyper-parallelogram in hyperspace formed by 7/ a 's. If 
d = 2, this region would be a parallelogram to the first region. Considering the Eq. (5.12) as 
a semi-linear one, we should use linear function and its minimal condition to distinguish the 
feasible region. The region is defined as: 

fc«(l - e- x ^) < e a < k a (l + e- x ^), 

in ref.[18] is a subset of feasible region in Eq. (5.12). According to appendix A, we have 
—k a < Pfs a < k a and by noting that J = NE and considering Eq.(3.4) we can write: 

(A - ]T P op )E = (EE Pafi)E a (5-13) 

P=0 a=l /3=0 

because E£ (a = 0, ...,d) are independent, we may conclude: 

Y,P op = N, e^ = ° 

/3=0 Q=l 

according (3.7), P a0 = 1 and by using relation (5.12) 

ka{ l _ e -x L m) < ^(k a - Pp a ) Va < k a {\ + e' x ^) (5.14) 

a 

which confirms the result of Ref. [18]. To investigate the abolition of asynchronization and 
also in order to calculate the scaling exponent of its suppression, we linearize the recursion 
map (3.11) near the fixed direction (synchronization state) of this map: 

lim n >00 Xi{n) = x (5.15) 

This leads us to write Xi(n) = x + 5xi(n). As expected, the maximum eigenvalue corresponds 
to synchronization state, giving the power scaling exponent of dynamical system[27]. Then, 
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by regarding Eq. (5.16) and by noting that Pslf is the second largest Floquite multiplier, the 
scaling exponent of suppression of the asynchronization, p is defined as: 

7J = In | ^ max |= — ln(second largest eigenvalue of A) = — In | \sle{A) | (5-17) 

PSLF 

where Xsle{A) means second largest eigenvalue of A. We can describe the transition proba- 
bility with respect to synchronized state by: 

II x N (n)-x syn || ~ e ~ nXsyn 

Transition probability changes in each time step (n — > n+1), and takes positive values if e > 0. 
The smaller the second largest eigenvalue modulus, the faster the Lyapunov exponent converges 
to its equilibrium distribution (Synchronization state). This problem can be formulated as a 
convex optimization problem, which can in turn be expressed as a semi-definite program [25]. 
Restoration of synchronization time tns is a time that in it system transits from synchronization 
state into asynchronization state. It is in fact the inverse of p 

tRS = ~ = ~ , , r 1 | - (5.18) 
p In | X max {A) | 

From (5.18) it is clear that, Synchronization time t is independent of Lyapunov exponent of 
map [29]. 



5.2 Stability analysis with Lyapunov function 

Lyapunov functional] is the generalization of an energy-like function that is decreased along 
trajectories. If the lyapunov function exists, then closed orbit will be forbidden. The forma- 
tion of synchronization can be studied through the Lyapunov function. It states that for all 
asymptotically stable systems, there exists a Lyapunov function whose derivative along the 
trajectories of the system takes a negative value. To study the stability region of the coupled 
map lattices, one can define lyapunov function in terms of E a as follows: 

Y J E a = I N -^ = E AS (5.19) 
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So, we can define the Lyapunov function as: 



N 



v ( n ) = J2 x i( n )( E As)ijXj(n) = x*(n)(I n -)x(n) 

i,j n 



(5.20) 



Lyapunov function takes " 0" value at synchronized state and in step n it becomes: 

N 

V» = $>(n) 2 



(E^^(n)) 2 



i=i 



N 



(5.21) 



Clearly V(n) > and the equality represents the exact synchronized state. For the asymp- 
totic global stability of the synchronized state, Lyapunov function must satisfy the following 
condition in the region of stability V(n + 1) < V(n) [32]. we can write of state in the following 
form x(n) — x s (n) + 5x(n) where 



x s (n) 



v 1 / 



x(n) 



(5.22) 



x s (n) is the synchronized state. By considering (4.9), V(n) = <5x*(n)5x(n), substituenting x(n) 
(corresponding to x(n + 1)) in (3.11) and by noting that 5x(n).x s (n) = 0, 



5x(n + 1) = $(x(ri)) eX + (1 - e)X' 8x(n) 



(5.23) 



according to (5.21), we can write the ratio of the Lyapunov function in step n+ 1 to Lyapunov 
function in step n 

2 

V(n + 1) 5^{n+l)5x{n + l) 6&(n) [eA + (1 - e)A"\ 5x(n)&(x(n)) 



(5.24) 



V(n) 5x t {n)5x{n) 
where according to reference [18], $(x) with N=even equal to By noting | Ai |> ... >| |, 
A = eA' + (1 — e)A" and according to appendix A, one can write 



5x t (n)A 2 Sx(n) \ < \\ Sx(n) \\ 2 max(\ \a | 2 ) 



therefore, 



V(n) 



< max{\ <t(x(n)) | 2 : x(n) E [0,1]}(| A \ 2 ) max 



< 1 



(5.25) 



(5.26) 
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The Eq. (5.26) is confirmed by the results of studies that have already been obtained in the 
literature [32]. The above mentioned result is converted to Lyapunov exponent method in 
section 5.1. For this purpose, we write the log of geometric mean of square root of Lyapunov 
function as: 



1 



n— 1 



lim — In i * ii I 

» W,, \ H'/') 



"tt 1 i V(m+1) ,\ ,. 1 TT - , 



n 



m=0 



\ 5f(0)5f(0) 



(5.27) 



by choosing 5x(0) as the idempotent of E a (a = 1, d), obviously, the righthand side of (5.33) 
equal to Al($) + In | \ a |, otherwise, if 6x^(0) E max Sx(0) ^ 0, it becomes + In | \ ma x |- 



Where AE„ 



XmaxEmax- It means that our results confirm Lyapunov exponent approach 



in subsection 5.1. 



5.2.1 Stability of periodic orbit 

In order to prove the stability of synchronized fix point we can rewrite the Lyapunov function 
and Lyapunov spectra at a synchronized state. By refereing to the definition of the coupled 
map lattice lyapunov spectra (5.7) at fix point x*, A^(x*) = In | \p& (x*) |. Both and 
| $'(x*) | are smaller than 1. Therefore, coupled map is stable. Also by refereing to the 
definition of the Lyapunov function at the synchronized state, we can prove the stability of 
fixed point at the synchronized state. We can rewrite relation (5.22) at fixed point as the 
following form: 



x s 



X" 



V 1 / 



(5.28) 



Where x* is coordinate of synchronized state. According to appendix A: 

V(n + 1) 



<| 4> 2 (x*(n)) |< 1 



V(n) 



(5.29) 
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by noting that x* is fixed point, then $(#*) is smaller than 1. So, for the stability of the 
coupled map the following condition should be met. 

V(n + 1) 



V(n) 



< 1 (5.30) 



5.3 Study of unstable synchronize and transverse state with Li- 
Yorke method 

Here we present the study of the topological effects on the dynamic behaviors of a coupled 
complex network in transition to chaos based on Li-Yorke idea. We considered a network 
consisting of nodes that are in non-chaotic states with parameters in non-chaotic regions before 
they are coupled together. It was shown that if these non-chaotic nodes are linked together 
through a suitable structural topology , positive Lyapunov exponents of the coupled network 
A/3 can be generated by choosing a certain uniform coupling strength e'(e"), and the threshold 
for this coupling strength is determined by the complexity of the network topology. 
Li-Yorke theorem [33] clarifies the existence of a snapback repeller, which implies the existence 
of chaos. This theory was generalized by Marotto to higher-dimensional discrete dynamical 
systems[34]. By employing Marottos theorem the transition to chaotic procedure is as follows: 
consider the N-dimensional difference equation 

x(n + 1) = $(x(ra)); x(n) E R n . (5.31) 

Suppose that the equation (5.31) has a fixed point x*. This fixed point x* is called a snap-back 
repeller if: 

a- $ is differentiate in a neighborhood B(x*,r) of x* with radius r > 0, such that all eigen- 
values of the Jacobian [D<&] are strictly larger than one in absolute values; 

b- There exists a point x(0) G B(x*, r), with x(0) ^ x*, for some integer m > 0, $ m (x(0)) = x* 
and $ m is differentiable at x(0) with det[D<5> m {x{0))] ^ 0. 



21 



If Eq. (5.28) has a snap-back repeller, then (5.31) is chaotic in the sense of Li-Yorke. At the 
synchronized state, the Lyapunov exponents of the N-dimensional dynamical system are 
described by (5.9). Now, by considering (5.9) at fixed point (x*) we have: 

Ap(x*) = In | \ p | +ln | $\x*) | (5.32) 

The transition to chaotic state takes place when: 

I A /3 \ W(x*y\ 

the second condition is: 

(m-l \ m m-l d 

1] &{x{i))) det(eA' + (1 - e)A" = J[ &(x(i)) J] A™ + (5-34) 
j=0 / i=0 /3=1 

Therefore, we need to choose coupling constants e' a and e' a so that all \p{f3 = 0, d) become 
nonzero. Also, x(0) is chosen in such a way that satisfies the $(x(0)) = and $'(x(0)) = 

or I 

As an example for the trigonometric chaotic maps [18], we introduce the above mentioned 
conditions: 

• For &n\x, a) (N = even), as it was discussed, [18] maps have only a fixed point attractor 
(x* = 0) then x — 1. 

• For $$(x, a ) (N = even), it was shown, [18] maps have only a fixed point attractor 
(x* = 1) then xo = 0. 

• By choosing odd value of N in '(x, a) , another fix point presented. For which the 
Li-yorke condition met too. 

The ratio of each polynomial volume to the to total defined volume is equal to "0" [35]. Then 
n«=i 7^ and Li-Yorke chaos is happened. 
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5.4 Kolmogorov- Sinai entropy 

Kolmogorov-Sinai (KS) entropy and Lyapunov characteristic exponents are two related ways 
of measuring disorder in a dynamic system. A definition of them can be found in many 
textbooks [20]. The KS-entropy of coupled map lattice can be assessed by computing Floquet 
multipliers integerals. It should be mentioned that computing the KS-entropy for measurable 
dynamical system ( such as (3.11) at synchronized state ) is also possible. The characterization 
of chaotic time dynamics can be made through the KS-entropy of the coupled map lattice. The 
ergodic choice of one-dimensional map x{n + 1) = $(x(n)) leads to the equality of KS-entropy 
with sum of positive Lyapunov exponents. Hence, according to Pessin's theorem [15], the 
ergodicity of one-dimensional map x(n + 1) = $(x(n)) implies the ergodicity of symmetric N- 
dimensional map (3.9) at unstable synchronized state ( synchronized state is stable for negative 
critical exponent A^). Obviously, the non-ergodic choice of x(n + 1) = <&(x(ra)) will lead to 
the non-ergodicity at synchronized state. The transition from chaotic synchronization (spatial 
order with temporal chaos) to non-synchronized states with positive KS-entropy (spatial and 
temporal disorder) occurs. Finally, comparing the KS-entropy with the sum of hp we have: 

KS($ N - coupled sychronized ) = ^ A /3 (<l> A r_ COU p, ed sychronized) . (5.35) 

k=i 

6 Some example of well known coupled map with cor- 
responding graph 

Globally coupled maps are one of the favorite models in the study of spatially extended dy- 
namical systems [36]. These models of the coupled map lattice correspond to the complete 
graph. A complete graph is a graph in which each pair of graph vertices is connected by an 
edge. The complete graph with N graph vertices is denoted and has N ^ N ~^ (the triangular 
numbers) undirected edges [38] . Globally coupled maps correspond to the trivial scheme with 
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agency matrices: 

A = I N , A 1 = I N - J N . (6.1) 

By considering Ei = I N — and E = and by using an arbitrary element of Bose Mesner 
algebra, the equation A = eA' + (1 — e)A" is reduced to: 

A = E + (1 - N^EJ (6.2) 

then eigenvalues of A are Ai = 1 and A 2 = 1 — N^- According to (5.18) and (A-2), it is clear 
that restoration of synchronization time tRs is reduced to: 

1 



tRS In | 1 - N% 



Based on (5.12): 



^ (1 - exp(-A L (x n+1 = $(x n )))) < 77! < ^(1 + exp(-A L (x n+1 = $(x n )))) (6.3) 

Similarly through Lyapunov function approach the following condition is met too: 

ki 1 ki 1 

|max$'(a; n )|) ~ Vl ~ | inax$'(a; n )|^ (6 ' 4) 

Coupling coefficients of the global coupled map (see as an example [37]) generates the agency 
matrix following (6.1,2). So, their stability analysis, stable region, can be calculated by re- 
lations (6.3,4). Also, in appendix C some of the well known graphs have been presented to 
simplify the generation of the new model. 



7 Conclusion and outlook 

In conclusion, it was found that the variety of the coupled maps (internal and external) display 
an identical pattern. Moreover, the synchronization capability of various kinds of coupling 
schemes were analyzed. We presented the proofs that the synchronization rate is mainly 
determined by the second largest eigenvalue of the Floquet multipliers. 
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In order to determine whether our results in this paper are generic or they only depend on the 
specific model, we also studied several other coupled maps. In future, the study may cover the 
following issues: 

• By applying the renormalization group theory in coupled dynamical systems, the scheme 
picture of model would be presented. As it was shown, near the synchronized state, 
different dynamical systems show the same patterns. This may lead us to this question 
"do the different dynamical systems follow the same pattern at synchronization state? 

• As it was shown in this study and in authors previous works [28, 29], in nonlinear 
maps, dynamical element represents the scheme picture, which may lead us to another 
interesting question " can we attribute the scheme to every nonlinear dynamical systems 
and coupled dynamical system too? " 

Further, we do hope that our obtained results through this paper will pave the way for further 
studies on nonlinear dynamical systems. 



Appendix A 

It was proved [39], that if A is a Hermitian matrix, then we have | x t Ax \ < max{\ Xa |) || x || 2 
and if x is the largest absolute eigenvalue then we have, max{\ Xa |), | x t A a x \ < k a \\ x || 2 , 
Thus for adjacency matrix A a , we have max{\ Xa u |) < k a and so Pp a < \ k a | for 
A = J2 a VaA a : 

I x l Ax |= I Va || x l A a x I < ^k a I r] a III x || 2 (A-l) 

a a 

Taking J2 a Vak a = 1 and by assuming that r] a is positive, we have: 

I x l Ax 



x II 2 



< 1 (A-2) 
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Appendix B 

We present some examples of the discrete one dimensional dynamical system with their invari- 
ant measure and Lyapunov exponent in order to simplify the generation coupled map model 
based on Bose-Mesner algebra: 



1. Bernuli shift map: 

<j>(x) = ^ (x - YllZl Pi) for E U Pi<x< Eti Pi 



(k=2,3,...) 



1 k 
fi=—, A = -Vpilnj9j 



2. Generalized Tent map 



\k+l k-l 



k-i 



= + (i + (_i)*) /or £ ft <x<'£ Pi (k = 2,3,...) 



i=i 



H=—, X = ~J2Pi ln Pi 
Pk tt 



3. Gauss map: 



= - - 



// 



A 



7T 



(l + x)ln2' 61og2 



4. Hut map: 



0(x) = ^ -1 + J9-16|x-|| 



^ = x + i A= I + 21n2-^ln3 



5. Piecewise parabolic map: 



(p(x) 



1 + r - yj(l -r) 2 + Ar \ l-2x 
2~r 
(1-r) 2 

H = l + r(l-2x), A = In 2 + ' ln(l - r) - 



/or r G [—1, 1] 

5r 2 - 2r + 1 , . , r 

In 1 + r + - + 1 

4r 2 
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Appendix C 

In this appendix we present some examples of the graph in order to simplify the generation of 
the coupled map model based on Bose-Mesner algebra. 

a- Strongly regular graph 

A regular graph G of degree r that is neither empty and nor complete is called strongly 
regular if every pair of adjacent vertices has exactly u common neighbors and every pair 
of non-adjacent vertices has exactly v common neighbors. The numbers r, u and v are 
the parameters of the graph. It is proved that a regular connected graph G of degree r is 
strongly regular if and only if it has exactly three distinct eigenvalues Ai = r > A2 > A3. 
If G is a strongly regular graph with parameters u and v then u = r + A 2 A 3 + A 2 + A 3 , 
and v = r + A 2 A 3 or 



A 2 , 3 = \[u-v± \j{u-v) 2 -±{v-r)\. 



then according to (5.18) and (A-2) it is clear that restoration of synchronization time 
tRS is: 

1 

tRS = : 1 = — 

In I I [u - v + yj (u - v) 2 - A(v - r)] | 

The parameters N, r, u and v determine not only each eigenvalue, but also the multiplici- 
ties of each eigenvalue. Since G is connected, the multiplicity of Ai is 1. The multiplicities 
m 2 and m 3 of A 2 and A 3 must sum to TV — 1 and by noting that the sum of the eigenvalues 
must be zero, the sum of the squares of the eigenvalues must be equal to iVAi = Nr 
.Thus m 2 and m 3 are 

\[y — 1 zp (A 2 + A 3 )(i/-l) + 2r 
2 A 2 — A3 

An example for strongly regular graph is Hamming graphs with diameter 2. H(2, q) is 
an SRG (g 2 , 2(q - 1), (q - 2), 2) with 



27 



Spec(H(2,q)) = 
if q = 2 then: 



Ai A2 A3 
1 mi m 2 J 



( 



\ 



2(9-1) 9-2 -2 
^ 1 2( ? -l) (9-I) 2 ; 



1 _ e --M$) 1 4. e --M*) 
: < 771 < : 



1 _ e -M*) 



< »7i + ^2 < 



1 + e- A ^*) 



2 -.*■«-- 2 

and if g > 2 using Eq. (5.14) we find a feasible region in the following form: 

1 _ e-M*) x + e~ AL(<I>) 
< 7/i + (q ~ 1)7/2 < , 

1 + (3 - 2q)e- x ^ 1 + (-3 + 2q)e~ x ^ 

9 — — 9 

then according to (5.18) and (1-4) it is clear that restoration of synchronization time t RS : 



t 



RS 



\n(q - 2) 



b- Cycle graph 

Another family of the N-coupled maps generated with cycle graph [38], is introduced as 
following: 

' 1 ... 

s= \ 1 ... 

... 

where S is the shift matrix. In cycle graph, adjancy matrices are introduced in the 
following form: 

N 

A = I N , A = S l + S-\ i=l, ...,[-], iV = odd 

N N 

A = I N , Ai = S l + S-\ i=l,..., 1, An=S~ N = even 

2 2 



V 



/ 
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Nearest neighbor coupled map lattice correspond to cycle graph. Nearest neighbor cou- 
pled map [36] is used as a model in the simulation of natural phenomena such as clouds, 
smoke, fire and water. It is one of the most important research areas in computer graph- 
ics. By using (C-6,7) and A = eA' + (1 - e)A": 

A=(l-Vi)I+j(S + S~ l ), 

where % — ee x + (1 — e)e' 1 '. Applying Fourier method we can write idempotents: 

i i N 



Ej 



where u> = exp(-rr), so eigenvalues are: 



Tj — 1 - 7/1 + rji cos 



2tu 
N ' 



The stability condition for cycle graph: 

I + e -Ai(X=*(x)) i _ e -\ L (X=*(x)) 

vi < on , ' for odd N > vi > for ever y N 

2(1 + cos f) 2(1 - cos f) 

According to (5.18) and (A-2) it is clear that restoration of synchronization time: 

1 



tRS 



In | cos ^ 
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